%test the admittance matrix
mpc = loadcase(case39);
Ybus = makeYbus(mpc.baseMVA, mpc.bus, mpc.branch);
Ybus = full(Ybus);

% 计算节点电纳矩阵中的shunt电纳矩阵
n = size(mpc.bus, 1); % 节点数
Ysh = zeros(n, n); % 初始化shunt电纳矩阵
for i = 1:n
    % 从case9系统数据中获取节点的电容和电感值
    C = mpc.bus(i, 9); % 电容值
    L = mpc.bus(i, 8); % 电感值
    
    % 计算节点的电纳值
    Ysh(i, i) = 1i * 2 * pi * 60 * C + 1i / (2 * pi * 60 * L);
end

% 打印结果
disp('Shunt电纳矩阵:');
disp(Ysh);